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1. Introduction: General structure of subtraction schemes 



In higher order calculations, the cancellation of infrared singularities is generally treated by 
the introduction of an infinitesimal regulator, eg in the form of a finite mass for the massless gauge 
boson, or by lowering the dimension D of the respective phase space integrals to D = 4 — 2 s. This 
way, the analytic cancellation of the respective divergent terms for fully inclusive variables after 
phase space integration is straightforward; however, numerical implementations of terms contain- 
ing small regulators prove to be challenging. In subtraction schemes, this problem is circumvented 
by the introduction of local counterterms, which mimic the behaviour of the squared real emission 
matrix elements in the singular regions; adding back the respective one particle integrated coun- 
terparts to the virtual contributions results in finite integrands for both real emission and virtual 
contribution phase space. Symbolically, this is given by 
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are the respective m, m + 1 phase space contributions to the total NLO cross section 1 . Convolution 
with jet functions then allows to define differential quanities and guarantees infrared safety of the 



respective Born contribution. In eq ( |1.2| ), the sum goes over all local counterterms needed to match 
the complete singularity structure of the real emission contribution. For each singular limit, the 
real emission matrix element factorizes according to 



|-^real \ 2 { P ) — ► D\^ Bom \ 2 (p), 



(1.3) 



where D denotes the dipole containing the respective singularity structure. As ^# rea i and ^#Bom 
live in different phase spaces, a mapping of the respective momenta from m + l to m particle 
phase space needs to be introduced, which is defined by a mapping function F map according to 
P = F mav {p). While the complete singularity structure of the limit considered is contained in D, 
both D and ~#Bom can depend on the leftover nonsingular parameters of phase space. D, and % are 
related by % = \i lE J dt, Dj, where the integration measure dt, is an effective one particle integral. 
Summarizing, any subtraction scheme needs to fulfill the following requirements: 

• definition of subtraction terms D, which, following eq (|L"3|), one by one mimic the behaviour 
of the real emission matrix element in each singular region such that their sum contains the 
complete singularity structure of the process, 

• definition of a mapping F map which guarantees total energy momentum conservation as well 
as onshellness of all external particles both before and after the mapping, 



'For hadronic initial states, an additional collinear counterterm rj c needs to be added in m particle phase space, 
which accounts for contributions already contained in the NLO PDFs. 
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• integration measure dB,, with a "smart" choice of variables providing optimal singularity 
structure parametrization. 

In the following, we will discuss this for our specific scheme, comparing with [jl|] where appropri- 
ate. 

2. Scheme setup 

In the scheme discussed in this report, the NLO subtraction terms are derived from the splitting 
functions introduced in [fj], and the m+1 torn phase space mappings needed correspond to the 
inverse of the respective shower m to m + 1 mappings. In the following, we will denote the m+1 
phase space four vectors by p\,p2, ■■■ and m phase space four vectors by p\, p2, .... Indices a, b 
will denote initial state particles; in m + 1 phase space, pi is the emitter, pj the emitted particle, and 
pic the spectator 2 . By default, for initial state emissions we set pi = p a in all general expressions 
below. Equally, we restrict our expressions to subtractions on the partem level and to the massless 
case; details on convolution with PDFs are given in and the extension to massive particles is in 
the line of future work. 

2.1 Momentum mapping and integration measure: Initial state 

For the initial state, the mapping from the m + 1 to the m particle phase space is given by 

Pa = [ 1 - P J ? I Pa, Pi = A(K,K)pi, p h = p b , (2.1) 

\ Pa-Pb J 

where the index / goes over all final state particles in the m particle phase space and with 

, — . 2(K + K)(K + K) 2KK 
A(K,K)=g v ; 1 ; + (2.2) 

where K = p a + Pb, Q = Pa + Pb, K = Q — pj. The phase space factorizes according to 

[d{p,f}m+l] g({P,f}m+l) = [d{p,f}m] d$ p g({p,f} m+ i), (2.3) 

where / denotes the flavour and with the D-dimensional integration measure 

2.2 Momentum mapping and integration measure: Final state 

For final state splittings, the initial state momenta remain unchanged; p a — p a , pb — Pb- 
The mapping uses all non-emitting particles as one spectator for momentum redistribution. We 
introduce the additional variables 

P = pi + pj, Q = pa+Pb, 

2 ra 



y = jp-^z-: a = ^f^I, A - x/n + v:^4,v. (2.5) 



In contrast to nip, in our case a spectator only needs to be specified if pj denotes a gluon. 
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The emitting particle is mapped according to 

_ P _ l-X+y 
^' X 2Xa 



Q. (2.6) 



All non-emitting final state particles are mapped using the Lorentz transformation as in eq. (|2.2|), 
where now K = Q — pi, K = Q — P. Especially, this means that the total number of mappings 
needed for a A^-jet final state scales as which reduces scaling 3 with respect to [|]] by a factor N. 
The phase space factorization takes a similar form as in the initial state splitting, ie we have again 

[d{pJ}m+\]g{{pJ}m+\) = [d{p,f}m] d£, p g({pJ} m+[ ) (2.7) 



where now 



x(2n) D 8w(p-X Pi - l -^ Q y (2.8) 

Jmax = (y/o — \fa — l) can directly be derived from total energy momentum conservation. 

2.3 Treatment of interference terms: dipole partitioning functions 

Double poles in splitting functions only arise if the emitted particle is a gluon; in this case, 
interference terms between different emitters have to be taken into account. In our scheme, we split 
the collinear and soft parts of the respective spin averaged splitting functions W according to [Q] 

Wu - W ik = (Wu - Wf) + (Wf - W ik ) , (2.9) 



where W u is the spin-averaged eikonal factor. The second part of the above equation can be then 



expressed in terms of dipole partitioning functions A' ik [|5j] 



Wf - W ik = 4 % a s A' ik rj^rro , (2-10) 

{Pj-PiPj-PkY 



where P& = pj ■ pipk — pj ■ pk Pi- Several choices for A' ik have been proposed [S]; all results given 
here have been obtained using eq (7.12) therein. 

3. Example of integrated splitting function: g — >■ qq final state splitting 

For our scheme, all collinear as well as singular parts of the soft splitting functions have been 
tested; a complete list will be given in [J6J] . In this section, we give the final state g — > qq dipole and 
the corresponding integrated term as an example, additionally commenting on the limit for m — > 2. 

For ag — ^ qq splitting in the massless case, the spin averaged subtraction term is given by D^|^#Bom(p)| 2 > 
with 



_ 4na s 

Dgqq — Tr — 



1 2 Z (1- Z ) 
1 -e 



(3.1) 



The subtraction scheme in pfl has a scaling similar to our scheme. 
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where we introduced the additional variables z = y=r , n 



l+v+A 



Q — jP, and all other variables as 
in eq. (2.5). This subtraction term was derived by squaring the respective final state shower splitting 
function in [0]. Momentum mapping is done according to eq. (E1J) and the Lorentz transform with 



the respective expressions for K, K. For the integrated splitting function, we rewrite the measure 



( |2.8| ) in terms of the variables introduced above and obtain 
(2 Pi Q) 1 - £ dQ. d -2 



which results in 



\6x 2 



(271 



l-2£ 



dzdyX l - 2e y- e 



[z ( 1 - z)] 6 ( Jmax - y) 9 [z (1 - z)] (3.2) 



2 16 2 
-3i"T + 3 [(a - 



1) ln(a — 1) — a In a] 



(3.3) 

As expected, for m = 2 the above expressions as well as the mapping completely reduce to the 
result in [jl]]. 

We want to comment that in our scheme, the most complicated expressions stem from the inte- 



gration of the interference terms as in eq. fl2.10| ). As all final state particles are mapped using the 
Lorentz transform A, the finite parts of the respective subtraction functions need to be evaluated 
numerically; details will be given in m\. 



4. First results 



As an example, we give the analytic result of our splitting functions when applied to dijet pro- 
duction at lepton colliders, as well as a numerical comparison for Drell-Yan process at NLO using 
[|l|] as well as the scheme proposed in this writeup. We additionally tested the scheme for Higgs 
production at hadron colliders and decay; the respective calculations will be presented elsewhere. 

4.1 Dijet production at lepton colliders 

For dijet production at lepton colliders, the final state squared splitting function D = D qqg is 
needed. We denote the four-momenta of the outgoing partons in this process with pi(q), Pi{q), 
and /53(g). The unintegrated dipole subtraction term for emission from p\ is then given by 
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with x n = 2 %® ■ The respective integrated averaged splitting function is 
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f) V Q 2 

Combining the above splitting functions for both emitters with the Born, real emission, and virtual 
matrix elements and integrating over phase space, we obtain the standard result 

4 



r NLO 



a NLO{2} + a NLO{3} = ^ 



2% 



-10+-7T 



/23 4 2 

+ 7T 

V 2 3 



T LO 



-^C F a L0 . 
4 n 



(4.3) 
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Figure 1: Relative difference between NLO corrections to single W production using Catani Seymour and 
Nagy Soper dipoles respectively, as a function of the hardonic cm energy. The results agree on sub-permil 
level, shown are the numerical integration errors. 

4.2 Drell-Yan production 

We calculated single W production for a pp initial state at NLO, using both the scheme in [jl|] 
as well as our scheme, including PDFs and varying the hadronic cm energy of the process. We 
here only show the numerical result for this process. Figure [| plots the relative difference between 
the two implemented schemes. We see that the numerical differences are on the permill level and 
consistent with zero. 

5. Summary and Outlook 

In this report, we introduce a subtraction scheme which reduces the number of mappings in 
the real emission part of an NLO calculations by a factor Afj ets with respect to the scheme 
suggested in We explained the setup as well as phase space mapping, and presented first 
results for an integrated splitting function as well as applications for simple processes. A complete 
listing of all integrated splitting functions as well as further examples will be given in 
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